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Evidence that a cluster of UHECRs was produced by a burst or flare 
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Abstract: The angular clustering of 5 Ultrahigh Energy Cosmic Rays (UHECRs) in the combined pub- 
lished AGASA-HiRes data has a probability of ~ 2 x lCF 3 of occurring by chance. A first analysis of the 
implications of the event energies and angular spreading is presented, which is applicable if the source is 
close enough that GZK losses can be ignored. Under this assumption, the observed energies of the events 
in this cluster favor a bursting rather than continuously emitting source, with the events emitted on a time 
scale short compared with 300 Dm P c years. Assuming the UHECRs experience many incoherent small 
magnetic deflections enroute from source to Earth, the arrival direction distribution allows estimation that 
{B 2 X)D ~ 7.7nG 2 Mpc 2 . If the spectrum at the source ~ E~ 2 , the total isotropic equivalent energy 
emitted in UHECRs is > 10 43 D^ pc ergs. 
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Introduction 

If UHECRs are charged particles and are produced 
by bursting sources, the events from a single burst 
should be clustered in energy as well as angle. This 
is because high energy CRs experience lower de- 
flection on average and therefore arrive to the ob- 
server more quickly on average than do lower en- 
ergy events. At any given place and time, the spec- 
trum has about a factor of 3 spread in energies. As 
time goes on, the average energy observed at any 
location decreases. Here, I examine a cluster of 5 
events in the combined published AGASA-HiRes 
data, whose spread in arrival directions is so small 
it could be consistent with instrumental resolution 
and is thus highly unlikely to be a chance cluster- 
ing. First, I compare the likelihood of fits assuming 
a continuous versus bursting source, then I assume 
the source is bursting and extract properties of the 
burst and intervening magnetic fields. The analysis 
presented here ignores energy losses during propa- 
gation; a more comprehensive analysis is in prepa- 
ration. 

The complete published data from AGAS A and the 
stereo HiRes consists of 57 AGASA events with 
nominal energy > 40 EeV and 271 HiRes events 
with nominal energy > 10 EeV. Both arrival direc- 



tions and energies of the former are published[5]; 
only arrival directions have been published for 
the latter[l]. With the HiRes collaboration, we 
searched for clusters of events in the combined 
high energy sample, 57 AGASA events above 40 
EeV and 40 HiRes events above 30 EeV, using a 
Maximum Likelihood technique [2]. One HiRes 
event at 38 EeV was found to be clustered with 
the known AGASA triplet, with an angular dis- 
persion consistent with measurement resolution. 
The probability of promoting the likelihood value 
of the triplet in the 57 event AGASA dataset, to 
that of the quadruplet in the combined 97 event 
dataset by chance, is 2 x 10 _3 [4]. This "promo- 
tion probability" is a useful indicator of the signifi- 
cance since it is not skewed by the existence of the 
original triplet, in case that had been a chance oc- 
currence. Other measures of the likelihood of the 
quadruplet being a chance association give a simi- 
lar significance[4]. 

A search of the remaining 23 1 HiRes events with 
energies in the range 10-30 EeV turns up a 5th 
event whose arrival direction is consistent with its 
coming from the same source with negligible mag- 
netic deflection or dispersion. Its association with 
the same source is less secure than for the original 
quadruplet, since there is a 1 in 6 chance of get- 
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ting as high or higher likelihood value as observed 
in the data when adding 23 1 events at random to 
the 97 original events. For the rest of the paper we 
proceed on the assumption that the 4 high energy 
events come from a single source, and will con- 
sider both cases that the 5th lower energy event is 
or is not from the same source. 



Energies favor a Bursting Source 

If UHECRs are produced by continuous sources 
and GZK distortions could be ignored, then the ob- 
served distribution of energies for a typical source 
would be the same as the distribution of all ob- 
served energies, and we could make a simple as- 
sessment of the probability that the cluster source 
is continuous, independent of knowledge of the 
exposure or the reliability of the energy measure- 
ments, since those affect the full dataset and the 
cluster equally on average. Namely, the probabil- 
ity of selecting - at random - 4 events above 30 
EeV and one or two events below, from a dataset 
with 97 events above 30 EeV and 231 events be- 
low. If the cluster is taken to include only the 4 
high energy events, the probability that their en- 
ergies are drawn at random from the full data is 
2 x 10~ 4 , while if all 5 events come from a single 
source the probability of finding the observed dis- 
tribution of energies is 10~ 3 . Due to GZK distor- 
tions this analysis can however be quite misleading 
as will be discussed in a forthcoming paper. 

In order to proceed in our analysis of the spec- 
trum of an individual source, we need to know the 
total exposure to the cluster. The stereo HiRes 
exposure is energy dependent and has not been 
published, although the stereo exposure has been 
shown in conference talks and the mono exposures 
have been published. The energy dependence of 
the stereo exposure should be very similar to that 
of a mono detector, so guided by figures shown at 
conferences, I normalize the integrated HiRes ex- 
posure, including weather cuts, to be equal at 70 
EeV to the integrated AGASA exposure, Sag — 
1500 km 2 yr, and I take the shape of the HiRes 
stereo exposure from the mono exposure using the 
functional form due to G. Hughes and D. Bergman 



(private communication), with the result 

£hr{E) = 7.510- 4 cxp[c(l-e(- a ( :E - b »)]km 2 yr, 

(1) 

where x = Log 10 (£;/10 16 eV). Above E = 
10 17 8 eV the parameters are a = 1.175, b = 0.732, 
and c = 14.893. While this is only an estimate 
of the HiRes stereo exposure, it is adequate for 
the purposes here. The total exposure is £ = 
£hr(E) + £ AG 0{E - 40EeV). The relative ex- 
posure as a function of position, rj, is normalized 
such that Jijdfl = A-k. In the direction of the 
cluster it is about the same for both AGASA and 
HiRes: r\ c = 0.2 sr" 1 . 

In a beautiful paper in 1978, Alcock and Hatchett 
obtained the distribution in arrival time and direc- 
tion of X-rays from an instantaneous point source 
and showed that it has a universal shape which de- 
pends on the source distance and parameters of the 
medium but not on details of the scattering process 
such as the form of the differential cross section[3]. 
Waxman and Miralda-Escude[6] (WM-E) adapted 
the Alcock-Hatchett analysis to cosmic rays from 
a bursting source undergoing multiple small mag- 
netic deflections. WM-E derive the flux of UHE- 
CRs with energies in the range {E, E + dE} re- 
ceived from a burst which produced N(E)dE cos- 
mic rays in the same energy range, by an observer 
at a distance D and time delay At (relative to pho- 
tons): 

3cE 2 N(E) 



(2) 

where q is the charge, A is the characteristic length 
scale of the turbulent magnetic fields and (B 2 \) is 
defined precisely in [6]. An expression for Gah, 
the normalized probability distribution function, 
and also the joint distribution in energy and angle, 
are given in [3]. Present statistics of this cluster 
are insufficient for the joint energy-angle distribu- 
tion to be useful. 

With an E~ 2 spectrum at the source, the shape of 
the observed spectrum (2) is shown in Fig. 1 as 
a function of E/E , ignoring energy loss during 
propagation. The parameter E is determined by 
the distance and time delay of the observation, and 
the magnetic structure of the intervening medium: 



E =D 



2g 2 (B 2 X) \ 
3cAt I 



1/2 



(3) 
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Figure 1: Observed spectrum for an E~ 2 spectrum 
at the source and an energy-independent exposure 
(solid) or the actual exposure (dashed), for sources 
close enough that energy losses during propagation 
are insignificant. 



The peak of the spectrum is at -E pca k = 0.214 E 
and the average energy is E — 0.276 E Q . 

The characteristic spread in arrival directions aver- 
aged over all events of energy E is 



a (E) 



^/2Dq 2 (B 2 X) 



3E 



(4) 



as for the case of a non-bursting source. Keep- 
ing in mind that the time delay At ~ D/c(l — 
cos(9 s ) 2 ) ~ D9 2 /(2c) , the form of the expres- 
sion (2) for the flux in each energy range can be 
easily understood: 



-fobs 



(Number of events produced at the so 
4ttD 2 (Spread in arrival times) 



si 



N(E) 



9cqN(E)E 2 



4irD 2 D6 2 /(2c) 4irD 3 Dq 2 {B 2 \}- 



and the apparent difference in normalization com- 
pensates the peak magnitude of Gah ■ 
Note that when observing an instantaneous source 
at a fixed time delay - as opposed to fixed energy - 
the distribution of deflection angles with respect to 
the direction of the source is flatter than a gaussian 
for small angles and drops faster at large angles. 
Waxman and Miralda-Escude[6] obtained an ap- 
proximate expression for the integral over energy 
of the joint AH distribution, at fixed time delay: 

p(se 2 ) ^ cM-0-58[se 2 /d s (E) 2 } 2 ), (5) 



where 59 is the angle between the arrival direction 
and the line of sight to the source. We will use 
equation (5) to fit for S (E). 

Given the exposure to the source and a spectrum 
at the source, say E~ p , we can find the value of 
i? which gives the best fit to the data as follows. 
Divide the energy range into N bins. For each bin 
the mean number of expected events is 

^(E ) = N(E ) E 2 -? G AH ((E/E ) 2 ) S(E), 

(6) 

where the normalization factor N(Eq) is chosen so 
that S^ 1 ^ i (i?o) = N c is the number of events in 
the cluster. The most probable value of E is the 
one which maximizes the likelihood measure 



(7) 



where V(/J,i, rij) is the Poisson probability of find- 
ing the observed number nj of events in the ith bin 
when fii are expected. We adopt henceforth p = 2 
for definiteness. This choice simplifies formulae 
and is likely to be close to reality, and the inferred 
E differs only slightly with other choices, e.g., 
p = 2.7. 

The energy of the 5th, low energy event has not 
been published, but HiRes has kindly released the 
event's ranking in energy (150th), constraining its 
energy to be within a few EeV of 15 EeV; this is 
plenty accurate for our needs. Using just the four 
highest energy events, the best fit value is E^ = 
190 EeV; with all five events E ( 5) = 150 EeV. 

The actual value of the likelihood measure defined 
above has no particular significance - for instance, 
it depends on the number of bins - but it can be 
used to assess the relative quality of different fits 
to the spectrum of the cluster by taking the ratio of 
the LM values. GZK distortions must be included 
for this to be meaningful for all source distances; 
the results will be presented elsewhere. 

Angular Distribution and Flux 

Henceforth we assume the events in the cluster are 
all protons. The best fit to the observed arrival di- 
rections under the bursting source hypothesis is ob- 
tained by maximizing the product of the angular 
probability densities for the N c events of the clus- 
ter, with respect to the direction of the source 6 
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and the energy-dependent magnetic smearing pa- 
rameter 9 S (E) introduced in (9). For each event, 
the 2-d probability density is proportional to 



exp 



-0.58 



exp 



{9i - 9 ) 



(8) 



This results in a source direction {RA, dec} = 
{169.67°, 56.70°} and 

6,(E) = 1.14°. (9) 

From these values one finds 

(B 2 X}D » 7.0 EcV 2 = 7.7nG 2 Mpc 2 , (10) 

where D is the distance to the source which pro- 
duced the quad. 

The normalization of the spectrum at the source 
can be estimated in units of {B 2 \)D 4 either 
by equating the total number of events, N c , or 
energy received from the source, E^-Ej, to 
/ E°' 1 F(E)£(E)dE, respectively. This gives 

N(E)E 2 3c 



= 2.4(...)10" 



(11) 



8ir(B 2 \)D 4 

respectively. Using the value of (B 2 X)D from the 
angular distribution gives 

E 2 JV(E) = 4.0 lO 42 ^ (j^j erg. (12) 

Integrating this over energy yields the isotropic 
equivalent energy the source emitted in ultrahigh 
energy cosmic rays with energies in the range 10- 
300 EeV: 

2 



E 



10-300 



= 1.310 



41! 



0s 

1.14° 



£>^ pc erg. (13) 



If the UHECRs from the source are beamed toward 
us in a cone of solid angle ACl, the total energy in 
UHECRs should be reduced by the factor |p. 

Since GZK losses have not been included, (13) can 
be considered a lower limit to the energy in UHE- 
CRs at the source. 

Minimum Flare Duration 



The source can be flaring rather than bursting and 
lead to the observed spectrum, as long as the 
source's peak luminosity is much greater than its 
normal or quiescent luminosity, for a timescale 
short compared to At. 

Conclusion 

The analysis presented here, ignoring GZK distor- 
tions, provides a first, naive look at the implica- 
tions of the Ursa Major event energies and angular 
separations for the nature of the source. It serves 
as a guide to the more complex analysis required 
when the source is at a large enough distance that 
GZK losses are important; the analysis and simu- 
lations for the more general case will be presented 
elsewhere. If the source is close enough for this 
analysis to be applicable, the spectrum favors a 
bursting or flaring rather than continuous source. 
In that case, the total isotropic equivalent energy 
emitted in CRs in the range 10 19 - 3 1 20 eV is 
> 10 43 ergs and the duration of the flare is 

<< 300 Z?m P c years. 
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From eqns 3 and 9 and the fit to 6 S we have that the 
arrival time delay is 

M= l 6 *^ ~ 300£, M P cyr. (14) 



